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Abstract 

It is shown that the spin is naturally introduced into classical mechanics if 
the latter is formulated as dynamics of the phase space density. It is shown 
that the uncertainty principle, as the amendment in this dynamics, restricts 
possible spins, and in particular equation for the particle with the spin h/2 
is derived. Also equation for the charge with this spin is derived when elec- 
tromagnetic field is included. In one example it is shown that the modulus of 
the spin changes with the gradient of the magnetic field. 

I. INTRODUCTION 

There are a number of parameters that uniquely specify a particle, and one of them is 
the spin . Yet there is a fundamental difference between, say the charge and the mass of 
a particle and its spin. The former are arbitrary parameters that enter dynamics equations, 
but the spin is the property of the dynamics equation, the Dirac equation p|. Therefore, in 
principle its explanation should be traced to the principles from which the Dirac equation 
is derived, and in this respect the spin appears not as a fundamental property. Despite this 
observation the explanation for the spin withstood various attempts, so that nowadays it is 
considered a property of the particle that does not have classical description. Nevertheless 
it is worth mentioning one attempt to explain it as a classical effect Q. For the electron to 
have the spin it is assumed to be a finite sized sphere (all experimental evidence shows the 
contrary) with certain distribution of the charge on its surface or the volume (what keeps 
the charge together is avoided as an answer) that rapidly rotates (to get the value h/2 for 
the spin superluminal equatorial velocity of the sphere must be assumed). The inability of 
the classical model (it is almost impossible to think of another classical model for the spin 
0]) to explain the spin resulted in the generally accepted view that in the quantum world 
there are phenomena that cannot be explained by our experience, and one of them is the 
spin Therefore one does not make models of the spin, one uses the rules by which it 

is mathematically described. One should say, though, that there are various semiclassical 
studies of the dynamics of the charges with the spin, which can be called empirical rules. In 
the classical equations of motion one incorporates the term that represents magnetic dipole 
|^-|TD[ (which is point-like and manifestation of the spin) interaction with the magnetic field, 
without questioning the source of this dipole. 

In this work it will be shown that the spin has natural description from the classical 
principles, but it has the feature that its value is arbitrary. However, by assuming the 



amendment that classical principles should be consistent with the uncertainty principle 
the spin is fixed and has the value h/2. This assumption in classical mechanics has been 



elaborated and tested in a number of cases ||TTHI4|, and therefore the idea behind it will 



only be briefiy described here. One starts from the observation that the deterministic view 
of classical mechanics, which is based on the idea that given initial position of the particle 
in the phase space its future and past is entirely determined if the forces on it are known, 
is fundamentally wrong. Even in principle one cannot imagine making experiment that 
does not incorporate inaccuracies in the measurement of the position of particles in the 
phase space, and this means that one only knows probability of their whereabouts. If one 
accepts this view then the particles in the phase space are not represented by points, but by 
phase space (probability) densities, and their dynamics is governed by the Liouville equation. 
Implementing uncertainty principle makes a restriction on the possible phase space densities, 
which has significant impact on its dynamics. As it was shown, this implementation results 
in deriving quantum dynamics equation, the Schroedinger equation IT^. Among the various 



examples treated in this way were diffraction and tunnelling effects |]T5[, for which it was 
shown to be described by solving classical equations of motion. It will be now shown how 
the spin is described, what is its classical description and what some of its properties are. 
Although combination of classical dynamics with the uncertainty principle was done in the 
relativistic case [16], when it was shown that Dirac equation is obtained, a number of points 
were missed in that derivation. The source of the spin is not discussed in details, and 
in particular why one needs relativistic theory to describe it. Also relativistic theory was 
described for a free particle only, and it is imperative to include electromagnetic (EM) field. 
In this work only non relativistic theory is described, and it is shown that the spin is not the 
relativistic phenomenon. Furthermore interaction with the EM field is described, and in one 
example it is shown that the spin does not have universal value, but its modulus changes 
with the strength of the field, more precisely it depends on the gradient of the magnetic 
component of the EM field. When the homogeneous magnetic field is restored (or when it 
goes to zero) the modulus of the spin gets back its value h/2. 

The discussion starts with a simple classical model that has analogy with the spin, in 
order to show that the concepts leading to its true description are not alien in classical 
mechanics. 



II. SIMPLE CLASSICAL MODEL OF THE SPIN 



The arguments against the possibility of classical description of the spin are very convinc- 
ing and yet one can show that this concept is naturally introduced into classical mechanics. 
In the traditional classical treatment a particle is represented as a point in the phase space 
whose dynamics is determined by two parameters: its initial position in this space and the 
force that acts on it. However, if the particle is represented as a probability density func- 
tion p{f,p,t) in the phase space, which takes into account all uncertainties connected with 
determining its initial conditions, then its dynamics is solved from the Liouville equation 

^ + ^-V.-p +F- Vp.p= (1) 
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The mass of particle is m , its momentum is p and the force on the particle is F. The 
simplest example to solve from this equation is for a particle that is not moving and it is 
placed in the origin of the coordinate system. This means that at t = the value of these 
parameters are defined as the averages 

J dh d^p fpo{f,p)^ J dh d^p ^ po {r,P) = 

and if the phase space density is a function of the module r and p then the (average) angular 
momentum of the particle is zero. However the (average) angular momentum squared is not 
zero, which is explicitly checked for the phase space density 

Po(r,p-) = -^e-^V'^^-^>^ (2) 
when this quantity has the value 

<L^>= J d\ d^p {fxpf po {f,p) = ^a'^y (3) 

Therefore particle at rest in the origin of the coordinate system manifests angular momen- 
tum, whose average value is zero, but not because its modulus is zero but because one sense 
of rotation in the phase space density is equally probable as the other. Indeed if the average 
over one sense of rotation, say the positive and parallel to the x-y plane, is calculated then 
the resulting (average) angular momentum is 

<L>= J d^r d^p fxppo{f,p) (4) 

d^r J dpr J dpe J dp^rp^sm.9 pQ{f,p) ^ —z 

where the components of the momentum with respect to the vector r were used. 

That the previous finding is not only of academic interest will be demonstrated on the 
example that is quite realistic. It is assumed that a charged particle is placed in a constant 
magnetic field, with the assumption that before t = the field is zero, and then its time 
evolution is followed. The problem has a straightforward classical solution except that one 
step is often overlooked. When one says that a particle with the charge e is placed in the 
magnetic field at t = then its equation of motion is 

mdlr^^x{vxA) (5) 

where the tentative assumption is that initial position tq and velocity are determined 
before the field is turned on. This assumption is fundamentally wrong because during the 
infinitely short time interval around t — Q the vector potential changed from its zero value 

— * 

to A, which means that it had a step-like time dependence. This means that the equation 
of motion for the charge is 

m dlr = ~dtA+ ^ x (v x = ~A 5{t) + ^ x (v x A^ 
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and the velocity in a very small time interval from t = — etot = + e becomes 



vo A{fo) + O{e) 

mc 



(6) 



It follows that in the magnetic field initial position for the charge is the same as when there 
was no field, but the initial velocity is given by or the momentum Pq = Po~^^ (^o)- This 
observation has important consequence on the time evolution of the phase space density. Its 
initial functional form, which is determined before magnetic field is turned on, is po{ro,po), 
but its time evolution is determined as if the initial momentum of the particle has the value 
Po = Po — The equation of motion (||) is therefore solved for the initial Pq, and 

solution for the position is r = f{fo, Pq, t) while for the momentum is P = g{fo, Po, t). The 
phase space density is therefore a function of the form p{f, P, t) and in terms of the initial 
phase space density it is given by 

p(f, P, t) = Po [/(f, P, -t),g{r, P, -t) + -^A (/(r, P, -t))] (7) 

Angular momentum of the particle is (4) but p is replaced by P. 

The simplest example is when magnetic field is uniform of the strength and it is 



directed along the z axes. The vector potential from which it is derived is A 



ho 



z X r, and 



solution of the equations of motion is (momentum Pq has only x and y components) 



^ cPo . / eh^t 

ro + — sm 

eriQ \ mc 



+ 



cPn X Z 



ehn 



1 — COS 



ehnt 



Pqz 
m 



P = Po COS 



ehot 



mc 



Po X z sin 



ehot 



mc 



mc 

+ Poz z = g{fo, Po, t) 



t Z = f{fo,Po,t) 



From these trajectories one calculates time evolution of the phase space density from (^, 
and by assuming that prior to turning the field on angular momentum of the particle is zero 
then immediately after that it is given by 



d\ £Po fo X Po Po fo, Po + -A (fo) 
L c 



eho 
' 2c 



(fro (fPo fo X 
rfVo d^Po {xl + yl) po{ro,Po 



Po--A{fo) 
c 



Po {ro, Pq 



Therefore the particle acquires angular momentum, and because it has origin in the internal 
rotation of the phase space density and not in the motion of its center of gravity, it can be 
shortly called spin. 

At any later time the angular momentum is 



L 



d\ d'P rxPpo fir, P, -t),g{f, P, -t) + -A f{r, P, -t) 
I c 



and after the change of integration variables its final expression is 



L = -z 



ehn 



Ac 



cos 



eho 



1 — cos 



ehot 
mc 
ehot 
mc 



j d\ d^P [x^ + y^) Po (r, P) - 
j d'r d'P [Pi + P,2) Po (r, P) 



(9) 
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It is a time dependent quantity, which is always oriented parallel to the z axes. 

The same example is calculated from quantum mechanics, in order to show that the 
result is not pure academic exercise. The assumption is that the initial conditions are the 
same, which means that the initial wave function /o is such that 



The function /o (r) is real because it represents particle that is not moving. Schroedinger 
equation is 



2m \ he 

and in the Cartesian coordinates it is (only rotationally symmetric problem is considered, 
in which case linear term in vector potential can be omitted) 



2m 



f 



In the z direction equation is for a free particle while in the x and y directions it is for har- 
monic oscillators. Both solutions are known, and the wave function is written as expansion 



f{x,y,z,t) = J2 J dkan^s{k)fn{x)fs{y)e'^"e 

n,s 



where fn{x) are eigenfunctions for the harmonic oscillator. The coefficients are determined 
from the initial /o, and because of its symmetry the summation indices run over only even 
integers. In the simplest model 

in which case the coefficients are 
where 

_ {2ehochf'^ fd^ {-2hc + ehodPf'^ 

On — 



(f)! V 2"-i(2/ic + eM')^"'''^^' 

From this solution one would want to calculate the (average) angular momentum, which 
is defined as L — m J d^r f x j where j is the probability current, and for a particle in the 
magnetic field it is given by 

J=-Im(rV/)- — 1/7 (10) 
m mc 
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It can be shown that the contribution of the first term in the current is zero, and so the 
angular momentum is 



Acehnd? 



Ah^c^ + e^/ioc?" + cos 



2!,2 .4 



ehot 



cm 



2„2 



(11) 



In the phase space density (|^) the parameters a and h are arbitrary but for legitimate 
comparison with quantum results they must be specified so that 



Po {r,p) = \f{r)f \g{p)\' 



-e 



(12) 



where glp) is the function /(r) in the momentum space. For this choice of the initial phase 
space density the angular momentum (P) has identical value as ([Tl|) , which is shown by 
evaluating simple integrals. Therefore classical and quantum dynamics of a charge in the 
homogeneous magnetic field produce identical solutions. 



III. THE SPIN 

There are two important results of the previous analysis for the following discussion. 
One is that if classical dynamics is formulated from the Liouville equation then its solution 
is identical to the quantum. It could be argued that the identity is accidental, but it will 
be shown that this is not the case. Second result indicates that there are two contributions 
to the total angular momentum of the particle. For a completely general initial phase space 
density it can be shown that the total angular momentum is given as a sum 

J = L + S (13) 

where the first term is angular momentum which is associated with translation of the centre 
of gravity of probability density P(r, t), and it is called the orbital angular momentum. The 
second contribution comes from the internal rotation of the probability density and can be 
called the spin. In the previous section it was shown that it is caused by the act of turning 
magnetic field on, because rotating electric field is induced. The spin is therefore not an 
unusual feature of the particle, but in that example its value depends on a number of external 
parameters, i.e. it is essentially arbitrary. One result in particular is intriguing, and it is 
summarized in (4). It was found that for a phase space density that represents stationary 
particle, but with spherically symmetric spatial and momentum probability densities, the 
average modulus of the angular momentum is not zero, but the angular momentum itself is. 
The question is can a phase space density be found with the same symmetry properties for 
which the angular momentum is not zero? Indeed it can be found, for example if one writes 

p{f,p) = po{r,p) + ^ (V^ X Vp) ■ s po{r,p) (14) 

then < L >= s but < > is not affected by this additional term. However, this phase 
space density is arbitrary because there is nothing that prevents choosing arbitrary value 
for the spin s . 
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The identity of the classical and quantum dynamics in the previous section is not acciden- 
tal, but a consequence of selecting phase space density that is in accord with the uncertainty 
principle. The principle states that if the phase space density produces standard deviation 
Ax for the Cartesian coordinate x, and standard deviation Ap^. for the momentum compo- 
nent Px then at each instant in time the inequality 

A,T /\px ^ c 

should be satisfied. Discussion how the constant c is determined from experiment is omitted, 
but anticipating value is c = h/2. The phase space density with this property is obtained 
by parameterisation (in the Fourier analysis it is called convolution, but it is also known as 
the Wigner function [p^^-pO|) 



p{r,p.t) 



1 



d'qe''^^-'^/^nf+q,t)f{f-q,t) 



(15) 



where f{r,t) should be determined by requiring that p{f,p.t) satisfies the Liouville equation 
(H). For a free particle this substitution results in the equation 



2ip-q/h 



(fqe 



which is satisfied if 



r(f + q,t) dtf{r- q,t) + ^nr+ g',t)A/(f - q,t) + 
f (f - g, t) dtrif + t)-^f{f- g, t)Ar (r + g, t) 



dJ{r,t) + -—Af{r,t)=0 
2tm 



(16) 



In the equation for the function f{f,t) one recognizes Schroedinger equation for a free 
particle. It can be shown, relatively easily for the harmonic type potentials, that the same 
identity holds for a general potential ||T^. Therefore the identity in the previous section 
is not accidental but result of imposing uncertainty principle on the phase space density. 
The only accident is that for the Gaussian type probabilities the parameterisation (^) is 



identical to (pIS]). 

There is one very important consequence when the uncertainty principle is imposed on 
the phase space density. If the (wave) function / is assumed to be spherically symmetric, 
i.e. it depends on the modulus of r, then it can be shown that no-matter what its functional 
form the value (^ for the angular momentum squared is always < >= |. Likewise it can 
be shown that although the angular momentum (4) is zero its value from only one half of the 
phase space is always < L >= (from now on h is set to unity). These are like universal 
numbers, in which case one can think of possibility to have phase space densities whose spin 
is not zero, always having certain universal value. To find the phase space density with that 
property one starts by calculating probability current from ([T5|) , and by writing momentum 



as 



p = S(S-p) +Sx (px S 



(17) 



it is given by 



— Im 

m 



rs{s-\/f 



— Im 

m 



riVxS) xSf 
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where S is some unit vector. Angular momentum that results from this current is zero, but 
if additional current is added that results from the second term in (JT^) one obtains 



j = — Im 

m 



rsis-vf 



Im 

m 



r{vxs)xsf + — {vxs)rf 



and by assuming that the last two cancel the resulting current produces angular momentum. 
This current must obey the continuity equation 



dt (/7) = - V ■ ^' = -- [/* ■ Vj (^5 ■ V j / 
and when equation (0) is taken into account then it can only be satisfied if 



For arbitrary S this identity is not possible, which means that from the simple parameteri- 
sation (|l^) one cannot obtain equation with the property that it exhibits the spin explicitly. 
The phase space density ([T5|) can be generalized as 



p{f,p.t) 



^ J d'qe''P^-^rif+ q,t)fif- q,t) + g*if+ q,t)gif- q,t)] 
^ / £q e^'P-'^F+ (f + g, t)F{f - g, t) 



(18) 



which is also in accordance with the uncertainty principle but allows much greater freedom 
to satisfy previous requirements. Important step in this goal is to note that if one writes 



F{r) = f{r) 



cos f e"'^<^/2 
sin I e*'^/^ 



/(f)f/(e,0) 



(19) 



where /(r) is a scalar, then in the product 



U+AU = cos^ ^ + ^2,2 sin^ ^ + (Ai,2e-'^ + ^2,16^'^) sin ^ cos ^ 

the matrix elements Ai^j can be chosen so that every Cartesian component of the unit radial 
vector f can be reproduced. It can be easily verified that with the matrices 



A = a^ 



1 

1 



A 



-i 

1 



A 



1 
-1 



any unit vector can be written as f = cr^^x + ay y + az z = a, where the product with the 
matrices and U is omitted but it is implied. The symbols cr^, ay and a^ are used because 
they traditionally represent Pauli spin matrices, but this is also known as the Cayley-Klein 
parameterisation ||2T| in the classical mechanics of rigid body dynamics which was developed 
much earlier that the concept of the spin was introduced into physics. By using spin matrices 
it can be shown that any vector, and in particular momentum, can be decomposed as 



P 



{a ■ p)a — i a X p = a{a ■ p) + i a x p 



(20) 
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where from now on the angles in the matrix U give direction of the spin. This is alternative 
decomposition of the momentum to that in ([TtD , which results in the probability current 



J = — Im [F+(f, t)a a ■ VF(f, t)] - -^V x F+(f, t)aF{f, t) (21) 

It can be shown that the angular momentum resulting from the current (^T]) is zero, but 
if the phase space density is modified as in (]T^), which in this case isQ 

p{f,p.t) = ^ J (fqe^''^-''F+{r + q,t)F{f~ q,t) + (22) 
(V, X Vp) • / Sqe^'^'^F^{f+ q,t)aF{f- q,t) 

27!-^ J 



where s is for the moment arbitrary, then the current (|21| ) is modified as 



^1 1 — s 

j = — Im [F+(r, t)a a ■ VF(f, t)] - — — V x F+(r, t)aF{f, t) (23) 



and the resulting angular momentum is not zero. The current ( PB| ) satisfies the continuity 
equation because the divergence of the second term is identically zero, while for the first 
term one can prove that (a ■ V) (c? ■ V) = A and F satisfies the equation ([T6|) . Therefore 
essential for the continuity equation is the first term in ( p3D but from the equation (|T6|) it 
is not possible to derive it. This is possible if the equation (|l^) is written in the equivalent 
form, as the set 

dtF = -a-VG ; G= -^a ■ VF (24) 

when it is simple to show that the first term in the current (p3D is derived. In this way the 
equations are derived that explicitly incorporate the spin of the particle. 

Still the parameter s in (p3|) is arbitrary, but its value is determined from the requirement 
that the total angular momentum of the particle should be conserved. From the time 
derivative 

dtL = ^dt j d^r d^pfxp j d^q e^'P'^F+{f+ q, t)F{r- q, t) 



where the function F satisfies the set of equations (|2J), and after somewhat lengthy but 
straightforward manipulations, which involves the set (|2^) and (^), one obtains 

d^L = -^dtj d^rd^p j d^qe^'P '^ F+{r + q,t)aF{f-q,t) 

which means that the total angular momentum can be defined either as 



^It should be pointed out that this addition to the phase space density in no way affects consistency 
with the uncertainty principle. 
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J = J d^rd^p j d^qe^'P'^F+{r + q,t) 



r X p -\ — cr 
^ 2 



F{r-q,t) 



or 



J = Jdrdpfxp p(f,p.t) 

where the phase space density is given by (^), but choosing s = |. This means that the spin 
of the particle is 1/2. In the first definition of the total angular momentum one integrates 
in all but spatial variables, in which case 



J = Re 



d^r F^{f, t) ( -f X V + ^CT ) F(f, t) 



from where one defines operator for the total angular momentum as 

J = -f X V + -a 
i 2 

This is recognized as the total angular momentum operator for the particle with a spin in 
the standard quantum treatment. 



It is important to make one observation. If one starts from the set of equations (p^ ) it is 
not simple to conclude that the spin of the particle is 1/2, because the continuity equation 
that one derives from it is 



dtF+(f, t)F(f, t) = V ■ Im \F+(f, t)a a ■ VF(f, t)l 

from where the only choice for the probability current is 



J = — Im \F+{f, t)a a ■ VF(f, t)] 

m 



There is no way one can deduce the remainder in (]23|), and in fact one would conclude that 
the spin is 1, if it is calculated as the average 

3 = 171 d^r f X j 



IV. CHARGE WITH A SPIN IN ELECTROMAGNETIC FIELD 



In the previous section it was shown how to obtain equation for a free particle which 
explicitly exhibits the spin. It is given by (12^) , but it is important to extend it for a charged 
particle in the electromagnetic field. One starts from the Liouville equation for a particle 
with a charge e 



dtP^ P- VrP + 

m 



W + —W{pA) - —ipW)A ■Vpp = 
mc mc 



(25) 
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where $ is scalar and A is vector potential, and W = —e V$ — ^dfA. In order to select the 
phase space density that is in accordance with the uncertainty principle, and that resulting 
equations are for the spin 1/2 particle, one starts from the parametrization 



p{r, p, t) 



d^q e2*?-(P+f ^^)F+(f + q, t)F{f- g, t) 



which is for a particle with no spin. It differs from (|1^) in one important aspect: the 
momentum is modified according to the recipe of Section II (see in particular equation 
(0) ), in order to accommodate its change when the charged particle is placed in the EM 
field. If this parameterisation is placed into the Liouville equation (p5D then by a set of 
transformations one finds the equation that F satisfies. They are quite tedious and lengthy, 
for example it can be shown that 



q- A 



VqF+{f+q,t)F{f-q,t) 



q- A 



F+{f+q,t)F{f-q,t) 



2i 

e 
c 



- / d^q e'^^^-^P+^^V, ■ V,F+(f + q, t)F{f- g, t) 



£q g2i9.(p+fA)^. V,F+(f + q,t)F{f- q,t) 



but one eventually ends with an expression for the Liouville equation that can be con- 
siderably simplified if it is assumed that the scalar potential $ is a quadratic function of 
coordinates and the vector potential A is linear. In such a case one can write 



q ■ V$(r) = -^{f+ q) - -<l>(f - q) 
and the Liouville equation is satisfied if F satisfies 



A{f) = ^l(f + q) + li(f - q) (26) 



i dtF{f, t) + 



2m 



V 



A F{f,t) +e^F{f,t) = 



(27) 



The equation for F is for the charge with no spin, and to include the spin one would 
follow the same procedure as in the previous section. However, this simple extension is not 
straightforward and before doing that it is necessary to discuss ordinary classical analysis 
of the time evolution of the phase space density ([I^) when EM force is included. 



A. Dynamics without uncertainty principle 

For a free particle with a spin the phase space density is given by (p!^, but when EM 
field is included then for t = its logical extension would be 

Po(r,p) = Pin{f,P+ + (Vr X Vp) Pin{rJ+ -A) (28) 
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and for arbitrary time 



p(f, p, t) = Pi. 



fJ+-A{f) 
c 



where the vectors / and g were defined in (H), but here they have functional dependence 
f{f,'p, —t) and glfjp, —t), respectively. The first term is not of immediate interest, however 
the spin term, the second one, is. It will be analysed for the simplest example, charge in 
a homogeneous magnetic field. The solution for the classical trajectories was obtained in 
Section II, summarized in (B), and therefore the spin can be calculated from 



S 



- I d-^r dyfx ps - (V/ X Vg)p^n 



f,g + -A{f) 

c 



To calculate the integral one needs to transform (V/ x V^) into the derivatives with respect 
to the coordinates and components of the momentum, which is straightforward from the 
knowledge of the vectors / and g, but the result is somewhat complicated. In short, one 
gets for the spin 

5" = [sa; + Sxy cos(a — ojt)] X + \sy + Sxy sm{a — ut)] y + Sz cos(co't) z 

where s^y is projection of the spin sin the x-y plane. Its typical feature is that its modulus 
is not preserved, but precesses in the x-y plane with the frequency u = ehQ/{mc), and 
its z component also oscillates. The solution appears acceptable, except that it does not 
conform with the notion of the spin should be, e.g. it is expected that its modulus should 
be preserved. Namely, if the spin should have any meaning then its value must not depend 
on the circumstances as described. It can be shown on another example that this definition 
of the spin suffers from further serious drawbacks, e.g. if perpendicular to the magnetic field 
one applies say a constant force then the time dependence of the spin, as it is defined, would 
have a more complicated form. Yet the spin, if it is defined as a property of the phase space 
density that depends on the magnetic field only, should not be affected by such translation. 
Therefore a more elaborate time dependence of the phase space density is required, which 
is more complicated than the one based on propagating 
One possible generalization of the phase space density is 



p{f,p,t) = po{r,p,t) + - (Vr X Vp) ■ [s{f,p,t)po{f,p,t)] 



(29) 



where s(r, p, t) will be called the spin function, and the spin is now defined as the average 



S 



d^r d^p s(r, p, t)pQ{f, p, t) 



This parameterisation restricts considerably the possible phase space densities, and its jus- 
tification is possible a posterior, i.e. when the uncertainty principle is introduced. Therefore 
in classical mechanics there is no justification for such parameterisation, but it is entirely 
consistent with the classical principles, as long as one finds the dynamics equations that 
produce it. The parameterisation (|29|) is possible at all times if the equation for the spin 
function is found, and to find it simpler problem will be considered: charge in the coordinate 
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independent magnetic field. In this case the spin function is only time dependent, and the 
total time derivative of the phase space density is 



1 e 
dtp{f,p,t) = -dts{t) ■ (V,. X Vp)poif,p,t) + - — s{t) ■ 

2 2mc 



Po{r,p,t) 



where it was taken into account that pQ{f,p,t) satisfies the Liouville equation 



dtPo H p- VrPo + 



m 



W H p X H 

mc 



VpPo = 



(30) 



and that dt (Vr x Vp) = Vr x [Vp x H j + (Vr- x Vp) dt (the proof of this is not given, but 
can be checked relatively easily). By slight rearrangement of the terms one gets that if 



dts(t) = — s(t) X H 
^ ^ mc ^ ^ 



(31) 



then 



dtp{f, p, t) 



2mc 



m 



Vp X ( X 



Po{r,p,t) 



which violates the basic principle on which dynamics of the phase space density is based on. 
i.e. the right side should be zero (Liouville condition). This, however, appears a difficulty 
but in fact it is not, if certain conditions are met. The basic one is that the total phase 
space density is conserved, which means that 

d^r d^p dtp{f,p,t) = 

This is indeed the case, but it is not the only condition. For example, from the integral 



d^r d^p p dtp{r, p, t) 



2mc 



d^r d^p ps{t)- VpX [H X Vr ) po(^, P, t) 



one would have to obtain the same equations of motion as the ones from which time de- 
pendence of po{r,p,t) is calculated. Detailed analysis shows that this is indeed the case. 
Likewise, for the angular momentum one would obtain 



d^r d^p fx p dtp{f,p,t) = dtsit) s(t) x H 



-s{t) X H 



2mc 2mc 

from which the spin equation (|3TD is obtained, which shows self consistency of the extended 
Liouville condition. 

From the previous analysis one obtains that dynamics of the phase space density is 
obtained from the set of two equations: one is (|30D and the other is (^TJ). The problem, 
however, is how the extension is made when the magnetic field is coordinate dependent? 
Analysis is more complicated than the previous one but the idea is the same, and the final 
result is the set of equations for dynamics of a particle with the spin 

dts{t,f) = -^s{t,f} X H (32) 



dtp= W + —px H+—Vh (s-H 

mc mc V 

from which the Liouville equation (^) modifies into 



mc 



dtpo H P- VrPo + 

m 



W + —px H + —Vh\s- H 

mc mc 



■ ^pPo = 



where the index of the gradient indicates that spatial derivatives are taken only for the 
magnetic field (the spin function is now coordinate dependent). 
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B. Dynamics with uncertainty principle 

Previous analysis defined dynamics for a particle with a spin in EM field, but important 
ingredient was omitted: it is not consistent with the uncertainty principle. In order to make 
it the phase space density should be parametrized as 

p{r,p.t) = j d^qe'^''^<P+-c'^^F+{f+q,t)F{f~q,t)+ (33) 

^ (V, xVp)- j d^qe^"^<P+-c^^F+{f+ q,t)aF{f- q,t) 

which is a straightforward generalization of (^21) . The function F{f,t) can be parametrized 
as (|T^, but now having more general functional character 

F{r,t) = f{f,t)U{f,t) (34) 

where U~^{f,t)U{f,t) = I = unit matrix. However, analysis simplifies considerably if again 
it is assumed that magnetic field is coordinate independent, which is also in accordance with 
the restriction (EBI). In this case one can write 



p{f,p.t) =1 I d'qe''^^<P^+'^''^r{f+q,t)f{f-q,t) + 

(V. X Vp) ■ m [ rf3ge2*^^-(^^+f^^)r(f + q,t)f{r- q,t) 



{a ■ H)a - a{a ■ H) 



7^3 
1 

which is in the form (pQl). The equation for the spin function is (pTD, which implies that 

dtU^it) aU{t) = —U+{t)a x HU{t) 
mc 

From (|20|) one can write 

a X # = — 
2i 

in which case the equation for the function U (t) is 

i dtUit) = —a ■ H U{t) 

2mc 

while that for the function fif^t) is (P7|). The product of the two functions satisfies the 
equation 

1 f ?p — \ p — * 

i dtF = V A] F H ■ aF + e$F(f, t) (35) 

2m \ c / 2cm 

It can be easily shown that this equation is derived from the set 



dtF = -a ■ (v - -A] G + e$F(r, t) ■ G = -^a ■ (V - -A] F (36) 
\ c / 2zm \ c / 

which is straightforward generalization of the set (0). 

Generalization of the equation ( ^5]) for the magnetic field that is not homogeneous is 
straightforward, but connection with the classical set of equations is lost, because of the 
condition ( pBf ) which is not satisfied. Therefore, the connection is only approximate, and to 
what degree will be checked on one example. 
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V. GAUGE INVARIANCE AND TWO EXAMPLES 



One very important aspect of the previous analysis was entirely neglected, and it is im- 
perative to mention it. The gauge invariance was not checked, which is a very important 
condition that should be satisfied. In the standard quantum treatment the problem ap- 
pears to have been solved, but simple example indicates the contrary. For example if one 
would want to obtain momentum distribution for a charge in, say, static magnetic field then 
according to standard definition 

9{p) = j dh f{r) e-'^^-P 

which is not gauge invariant. One could try various amendments to obtain the invariance 
but non of them are satisfying. In this work gauge invariance comes naturally, and the 
problem of the momentum distribution is relatively easily solved. One should go back to the 
observation in Section II that momentum in the EM field is different from its value before 
the field is turned on. The relationship is summarized in (j^) (for simplicity static EM field is 
considered), which means that because momentum p is gauge invariant then momentum P 
is not. It follows from the relationship that in order to satisfy this requirement then the 
gauge A A + VO implies P P — VO. This, however, implies that the parametrization 
(PSP is gauge invariant (small p is used for convenience, but it is in fact capital P), and that 
under the gauge transformation the equation (pSl) transforms as 



1 / • • • \ 2 



idtF = V + - Ve A Ve F H ■ oF + e$F(f, t) 

2m \ c c c J 2cm 

and it is invariant. However, instead of transforming momentum one can formally assume 
that the wave function transforms as F — > Fe^®, in which case gauge invariance is satisfied 
without assuming transformation properties of V. 

One example in which momentum distribution will be calculated is generalization of the 
problem that was treated in Section II. In addition to charge being in homogeneous magnetic 
field it will be assumed that it is also in the field of three dimensional harmonic oscillator, 
and before switching on magnetic field the charge is assumed to be in the ground state of 
that oscillator. The problem will be solved by classical equations of motion, and from the 
previous analysis it follows that the result is identical to that obtained by solving equation 
(E^). In this example dynamics of the spin is not analyzed, only the probability densities 



The equations from which dynamics of the phase space density is deduced are 

dtp= -\7^ + —px H + —Vh (s- h) ; dtsit) = —sit) X H 
mc mc V / mc 

and for homogeneous magnetic field the last term in the first equation is zero. The potential 
for harmonic oscillator is. $ = ^muo^r'^, and if one does re-scaling s/h ^ s, ctx t and 
rx — i> r, where x = mc/h, then the first equation is (the spin equation is not analyzed) 
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d^r = ijj dt? z — uJq r 



where now oscillator frequency ojq stands for uiqIiI {mc^) and the cyclotron frequency uj for 
eHoh/ {m?c^). For simplicity this equation is solved in the x-y plane (magnetic field does 
not have effect on the motion along the z direction), and the solution is 



r = tq 



cos ( y 1 COS 1-^/^2 + 4^0^ ) + 



tuj 



+ 



Z X Tq 

2vo 



tuJ 



sin — cos -J^^ + Au^ - 



UJ 



too 



COS — sin -Juj^ + Au^ 



+ 



cos I ^ I sin 1 X A /coi^ + 4a;o 1 + 



2zXVq ■ ft^\ ■ f t / 2 , a 2 



where tq and Vq are initial position and velocity of the charge, in the x-y plane. According 
to prescription (J^) one obtains the phase space density at any time, and result is quite 
complicated. However, probability densities P(r, t) and Q{p,t) are simpler, and given by 



P(f,t) 



2 4 + g2 



exp 



2(4 + g2) ■ 

d^Ar 



4d^A + q' 
Q{P, t) = — X — exp 



TT 



4rf2(4 + g2) ■ 



where 



= 8 



+ q +q cos I txl^jT- 



A„ = 16 + lOg^ + - 2q^ cos t 



4: + q^ 



d4 



where q = u/ujq and d = l/y/uJo is the width of the ground state of the harmonic oscillator. 
Interesting limit is g ^ 1, in which case the width of 



P{r,t) 



1 



1 + cos (tuj) 



exp 



2r2 



c/2 [1 + COS (too)] 



oscillates between d (unperturbed oscillator) and (nearly) zero. On the other hand, momen- 
tum distribution has the limiting form 



Q{p,t) 



exp 



4rf2 



-P 



which is (nearly) time independent, but it is very wide. 

In the second example dynamics in the inhomogeneous magnetic field of the charge with 
a spin will be analyzed. For the vector potential it is assumed 

A = hiz {—y X + X y) 



and the magnetic field is 
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— * 

H = ~hix X — hin y + 2hiZ z 



In the scaling defined earlier the classical set of equations is 

; dts{t, r)=E s{t, r)x H (37) 



r= e 



while the equation (BSf) (quantum equation) is 



id^F = -ieAf F -^H ■ aF (38) 

where e = eh^hi/lm^c^), and in vector potential the coefficient hi is omitted. The equations 
(P7|) and (PBD, together with the appropriate initial conditions, determine dynamics of the 
charge. As it was shown, if the magnetic field is coordinate independent, but with arbitrary 
time dependence, the two sets produce identical results. However, it is not clear that the 
two dynamics produce the same result if this condition is not satisfied, and the purpose of 
the following analysis so to check if this were the case. 

In classical dynamics (which is based on the set of equations ( pT]) ) one calculates time 
dependence of the phase space density (|33|) , and from parametrization (|3^ ) one assumes its 
general form 

p{f,p,t) = po{f,p,t) + ^ (Vr X Vp) • s{t,r)po{f,p,t) 



where 

Po{r,P, 



t) = i- /" £qe^'^<P^'^^f*(f+q,t)f(r-q,t) 
■n^ J 



The dynamics of the particle is therefore determined if initial po{r,p) is specified together 
with the orientation of the spin s{t = 0, r) = sq (its modulus is fixed). From the phase space 
density various quantities can be calculated, but one in particular is the spin 



s{t) = J d r d p f X p Psp{r,p, ^) ~ J ^ P ^(^' Po{^iPi t) (39) 

In practical implementation one generates random initial r and p from the distribution 
\po{r,p)\ and solves equations of motion (^). For N randomly chosen initial conditions the 
spin, for example, is approximately (for N ^ oo the result is exactly the same as (^) ) 

^^^^^ _ En=i^(^>^") Sign[po{fn,pn)] 

where Sign [po(^n,Pn)] is the sign of the initial phase space density for the n — th pair of 
initial conditions. 

Solution based on the equation (|38| ) (quantum solution) is more difficult to obtain. First, 
because of the cylindrical symmetry the equation is transformed into the cylindrical coordi- 
nates, and if F is parametrized as 
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1 
e'^ 



Gn{r,z,t) 



n 



0,±1,±2,.. 



where r is the radial distance in the x-y plane and is the azimuthal angle, then the equation 
for Gn{r, z) is angle independent. If the initial F is assumed to have the form 



F^fo{r,z) 



sin ^ e^l'^ 



cos ^ e 



where a and /3 arc directional angles of the spin, polar and azimuthal angles, respectively, 
then it is obvious that single is not sufficient to describe dynamics. One needs Fq and 
F_i to describe it, and the initial values for Gnij'^ z, t) are 



Go(r, 2, 0) = v^/o(r, z) 



cos I e" 







G'_i(r,z,0) = v^/o(r,z) 







sin f e'^/'^ 



and the overall solution is F = Fo(r, 2;, t) + -F-i(r, t). In order to simplify analysis it will 
be assumed that the spin is initially oriented along the z direction, in which case only one 
equation needs to be solved, for the index n = 0. The equation for Go(^) t) is 



i dtGo = - 



1 



Go + ^ra^Go - ^ez [a^ + /) Go 



and it was solved numerically by replacing partial spatial derivatives with their finite differ- 
ence approximate forms (care should be taken in implementing this procedure for the radial 
coordinate, but the details are omitted). In this case the equation is replaced by a matrix 
equation 



i dtGo{rn, Zm, t) = ^0{r ■,rj,Zk) GQ{rj,Zk,t) 

3,k 



(40) 



which is solved numerically as a set of differential equations. The accuracy of the solution 
depends on two factors. One is the step that is used for approximating the derivatives, and 
the other is how well the boundary conditions are satisfied. The limitation on the former 
is the amount of computer time one has at disposal for calculation. Limitations from the 
latter are in the bounds that are set on the range of coordinates, because one chooses Zmin 
and Zmax for and rmax for r, within which the solution is confined at all times. Implication 
is that at these ends Go(^max, -2min,max, ^) = 0, which means that the solution is accurate for 
all times when this condition is satisfied. 
Explicit example was calculated for 

1 .2 

/o(r) 

and the initial phase space density is 

Po(r,p) 



e 
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Initial orientation of the spin was assumed along the z axes (a = 0). Numerical values for 
the parameters are li = 10 and e = 0.01, which were chosen so that the important effects are 
easily noticed, however they are not necessarily easily realized in experiment. Thus for the 
electron d ~ 2.4 * lO^^^m and hi ^ lO^^T/m. This field gradient is enormous and beyond 
means of any practical set up, but the effects of it on dynamics of the charge with the spin 
are clearly manifested. 

Two quantities were calculated from classical and quantum equations, one is the spin, 
which in classical dynamics is defined as (pQ]) and in quantum as 



m = ljd'r G+aGo (41) 

and the other is the probability density for finding the charge along the z coordinate, and it 
is defined as 

P{z,t) = J d^p dx dy po{r,'p,t) = J dx dy GqGq 

In the classical calculation 10^ trajectories were calculated to obtain reasonable accurate 
phase space density. In the quantum calculations the following parameters were chosen: 
Zmin = —100 and Zmax = 200 with the number of divisions of this interval riz = 600, and 
'"max = 200 with the number of divisions = 200 (total of 240000 equations in the set 
(^) ). Large value of these parameters to obtain (reasonable) accurate quantum results (as 
the rule for longer time evolution one needs considerably larger set of equations (^0]) ) is 
an indication that dynamics is dominated by some irregular process, about which it will be 
learned shortly. Time dependence of the spin, which is always pointing in the z direction, 
is shown in Figure 1, where the solid line shows quantum result and broken line classical. 
The agreement is nearly perfect up to the time t = 30 (~ 3.9 10~^° sec) when deviation sets 
in. The possible origin of this deviation will be discussed shortly, but the result in itself is 
interesting to emphasize. Both the classical and quantum dynamics indicate (in more or less 
identical manner) that the spin does not have a fixed value, but its modulus changes with 
interaction. In fact it changes with the gradient of the field, because as soon as the field 
becomes uniform the spin restores its original value 1/2. It is interesting, however, to note 
that the classical equation of motion for the spin function in { ^7\ j conserves its modulus, 
and the same is true for the quantum equation, but the spin, being either classical (|39|) or 
quantum ( ^I]) average over the spin function, it is not conserved. 

The probability density P{z, t) for the time interval till t = 20 is shown in Figure 2, where 
the solid line is quantum result and broken is classical. The difference between the two results 
is negligible. However, for the time beyond t = 30 it becomes noticeable, as shown in Figure 
3. The essential feature of the probability density is that it consists of two parts, the central 
and stable peak around z = and its extension in the positive z direction. The central peak 
is dominant, and the two calculations give nearly identical results, however, the extension in 
quantum calculations is (randomly) oscillatory while classical result appears as the average 
of it. Some of the oscillatory structure in quantum result is manifestation of numerically not 
converged calculation, but overall structure is reliable. From the analysis of the previous 
section it is expected that classical and quantum results are not identical, but the surprising 
finding is that quantum is so difficult to obtain and it has such a (nearly) random structure. 
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Explanation for the discrepancy in Figure 1 can now be given. The spin is the average over 
the spin function, whether classical or quantum, and because quantum probability shows the 
oscillatory structure and classical do not the averages cannot be expected to be the same. 
In fact the deviation between the classical and quantum calculation of the spin starts to 
be noticeable at the time when the quantum oscillations start to dominate the probability 
density. In any case the exact agreement is not expected, as derivations of the previous 
section indicates, and when this is taken into account discrepancy between classical and 
quantum results is not fatal, in fact it is quite good. So much so that it can be claimed that 
the spin has classical and non relativistic origin, although precise dynamics of the particle 
with a spin should be described by the set of equations (|36|). However, classical description 
of the spin has one great advantage, it offers understanding of the underlying processes that 
contribute to the final results, shown in Figures 1, 2 and 3. They are result of averaging 
over large number of trajectories, each with different initial conditions for the coordinates 
and momentum (the spin is always fixed). Few typical ones are shown in Figure 4, where 
time dependence of the z coordinate of the particle, the z-th component of the spin function 
and projection of the spin function in the x-y plane are shown. The modulus of the spin 
function is always constant for each trajectory, but its x-y projection is not zero. On the 
other hand, both quantum and classical calculations indicate that the spin has only the z-th 
component, which means that the zero x-y component of the spin is result of averaging. In 
fact in classical calculations with a finite number of trajectories the x-y component of the 
spin is not zero, but as its number is increased it converges to zero. Another important 
result is shown in Figure 4. Each trajectory is different, with no relationship between them. 
Even two very close with the initial conditions soon become unrelated, and indication of 
the underlying chaos in dynamics. This means that the "chaotic" oscillatory structure in 
Figure 3 is manifestation of the chaotic underlying dynamics, and explanation why there 
was difficulty to obtain reasonable accurate results. 



VI. SUMMARY 

Analysis of the spin, and calculation of examples, indicate that the spin is classical in 
origin, although in summary it could be labeled as the "fundamental parameter" of the 
particle, rather than derivable from some classical model like a spinning sphere. The reason 
why the spin should be regarded as such is in the way how the spin function was introduced. 
By parametrizing the phase space density as (PW), and enforcing that it retains this a form 



in the course of time, has as the result that the spin function is determined. It is a solution 
of an additional equation for the dynamics of the particle, so that a complete classical set of 
equations is (|3^ . The first equation is, without the spin function term, the usual classical 
equation of motion, and by including this term it is coupled to the second equation for 
the spin function. There is a very profound impact of this dynamics on the phase space 
density. The spin function is an additional degree of freedom for the classical particle, which 
means that at a given point in time two, or more, classical trajectories can cross (this is 
not possible for a spinless particle) in the phase space. This is possible because trajectories 
are not determined by only the position and momentum of the particle, but also by their 
spin function. Therefore if for a given time t two trajectories cross in the phase space this 
means that their spin functions are different. As the result the requirement that the value 
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of the phase space density is preserved in time, i.e. dtp{f,p, t) = 0, from which the Liouville 
equation for the spinless particle is derived, is no longer true. For the particle with the spin 
the Liouville equation is derived from the inhomogeneous relationship (i(p(r, p, t) = f{f,p, t), 
but in order to preserve the physicality of the problem the function f{f,p,t) must satisfy 
certain requirements. Those were discussed in the section following the equation (pQ]). In 
this respect the spin appears as a fundamental property of the particle, in addition to its 
charge and mass, even in classical mechanics. However, the observed value of the spin is the 
average over the spin function, and from this point it is a quantity whose magnitude greatly 
depends on the forces that act on the particle. 

There is another point that needs a brief discussion. It is customary to think not only 
that the spin does not have classical description, but also that it is inherently a relativistic 
phenomenon, i.e. result of making quantum theory relativistic, but in a specific manner as 
Dirac did. The original intention of Dirac was not to look for the equation that describes the 
particle with the spin. The intention was to make the Klein-Gordan equation first order in 
the space-time derivatives, for the reasons that can be found explained in any textbook on 
relativistic quantum theory. In doing so four component wave function was introduced, and 
the equations were found for their components. In the non relativistic limit these equations 
reduce to (^) for a free particle and to ( |5BD for a particle that is in the EM field. However, 
in this work the pursuit was to find equations for the particle with the spin, and it was 



explicitly shown that they can be obtained if and only if one makes generalization ([TsD of 



the phase space density. This is equivalent to requiring that the wave function is not single 
but the four component object that produces the set (|2^) , which is also first order in space- 
time derivatives, the same as the Dirac equation. In conclusion, therefore, one can say that 
Dirac derivation was accidental discovery of the spin, and because no alternative arguments 
(Pauli spin equation is an ad-hoc derivation, and can only be traced to the non relativistic 
limit of the Dirac equation) for deriving the Dirac equation were found the only realistic 
conclusion was that the spin is of the relativistic origin, with no true explanation in non 
relativistic quantum theory (not to mention classical theory). 
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FIGURES 



FIG. 1. Quantum (solid line) and classical (broken line) time dependence of the spin of the 
spin 1/2 particle in inhomogeneous magnetic field, whose axes of symmetry is z axes. Initially the 
spin of the particle is parallel to the z axes, and in the course of time it remains so. 

FIG. 2. Quantum (solid line) and classical (broken line) time evolution for the probability 
density P(z,t), up to the time when quantum and classical calculations for the spin of the particle 
start to separate, as shown in Figure 1. 

FIG. 3. Quantum (solid line) and classical (broken line) time evolution for the probability 
density P(z,t), beyond the time when quantum and classical calculations for the spin of the particle 
produce nearly identical results. 

FIG. 4. Typical time dependence of the z-th coordinate of the particle, the z-th projection of 
its spin and the projection of its spin in the x-y plane. Small changes in the initial conditions result 
in entirely different time dependence of these parameters. 
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